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Abstract 

The  Manned  Right  Simulator  at  the  Naval  Air  Warfare 
Center  in  Patuxent  River,  MD  maintains  high  fidelity 
fixed  and  rotary  wing  simulation  models.  The  aircraft 
simulations  are  utilized  for  a  wide  range  of  activities 
including  flight  test  support,  pilot  training,  and  control 
law  analysis  and  design.  Validating  airaaft  math 
models  against  flight  test  data  is  an  important  pan  of 
the  simulation  process.  Linear  model  comparison  was 
used  to  validate  the  lateral-directional  dynamic  modes  of 
the  V-22  tilt-rotor  aircraft  in  airplane  mode.  In  this 
technique,  linear  model  approximations  of  the 
simulation  and  aircraft  dynamics  are  calculated 
independently  and  then  compared.  The  simulation  linear 
state-space  model  was  extracted  from  the  nonlinear  V-22 
simulation  using  a  perturbation  method.  The  aircraft 
linear  state-space  model  was  fit  to  flight  test  data  from 
lateral-directional  maneuvers  using  parameter 
identification  tools.  Time  history  comparisons  were 
used  to  verify  both  linear  models.  Comparisons  of  the 
lateral-directional  modes  and  the  stability  and  control 
derivatives  of  the  two  models  were  made.  The 
differences  between  the  two  models  were  used  to  locate 
potential  problems  with  the  nonlinear  simulation. 
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Nomenclature 

aiL  Left  rotor  longimdinal  flapping  (deg ) 
ajR  Right  rotor  longitudinal  flapping  (deg) 
bjL  Left  rotor  lateral  flapping  (deg) 

biR  Right  rotor  lateral  flapping  (deg) 

A  State-space  system  matrix 

B  State-space  control  matirix 

C  State-space  output  mauix 

C]  Nondimensional  roll  moment  coefficient 

Cn  Nondimensional  yaw  moment  coefficient 

Cy  Nondimensional  side  force  coefficient 

D  State-space  output  control  matrix 

g  Gravity  constant  (ft/sec^) 

G  Cost  function  weighting  matrix 

Ixx  Moment  of  inertia  about  x  body  axes  (slug-ft*) 

Iyy  Moment  of  inertia  about  y  body  axes  (slug-ft^) 

Iz2  Moment  of  inertia  about  z  body  axes  (slug-ft^) 

Ixz  Product  of  inertia  about  x-z  body  axes  (.slug-ft") 

J  ,  Cost  function  matrix 

L  Roll  moment,  body  axes  (ft-lb) 

m  Aircraft  mass  (slug) 

N  Yaw  moment,  body  axes  (ft-lb) 

p  Roll  rate  perturbation,  body  axes  (rad/sec) 

r  Yaw  rate  perturbation,  body  axes  (rad/sec) 

t  Time  (sec) 

T  Total  number  of  time  steps 

Tr  Roll  mode  time  constant  (sec) 

Ts  Spiral  mode  time  constant  (sec) 

u  State-space  input  vector 

U  X  velocity,  body  axes  (ft/sec) 

Undi  Theil  coefficient 

V  y  velocity  perturbation,  body  axes  (ft/sec) 

W  z  velocity,  body  axes  (ft/sec) 

x  State-space  state  vector 

y  State-space  output  vector 

Y  Side  force,  body  axes  (lb) 
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P  Angle  of  sideslip  (deg) 

63  Aileron  position  (deg) 

8,  Rudder  position  (deg) 

A  Denotes  perturbed  value 

Damping  ratio 

((>  Roll  attitude  perturbation  (rad) 

0  Pitch  altitude  (rad) 

p  Air  density  (slug/tf ) 

(On  Undamped  natural  frequency  (rad/sec) 

Subscripts 

0  At  reference  flight  condition 

Superscripts 

Time  rate  of  change 
^  Estimated 

Introduction 

The  V-22  tilt-rotor  simulation  at  the  Manned  Right 
Simulator  (MFS)  is  a  high  fidelity  nonlinear  simulation 
utilized  for  flight  test  support,  pilot  training,  and 
control  law  analysis.  The  airframe  model  is  based  on 
the  Bell  Helicopter  Generic  Tilt-Rotor  simulation.  The 
rotor  model  is  a  disk  model  with  dynamic  longitudinal 
and  lateral  flapping  states.  The  rotor  model  structure  is 
based  on  helicopter  aerodynamic  theory,  flight  test  data, 
rotor  test  stand  data,  blade  element  model  data,  and 
airplane  mode  propeller  efficiency  data.  The 
aerodynamic  model  is  a  component  buildup  of  the 
aerodynamic  subsystems  including  the  fuselage,  wing- 
pylon,  horizontal  tail,  vertical  tail,  and  landing  gear. 
The  aerodynamic  model  structure  and  data  tables  are 
based  on  wind  tunnel  data,  theoretical  equations,  and 
flight  test  data. 

The  V-22  simulation  i$  implemented  in  the  MFS 
Controls  Analysis  and  Simulation  Test  Loop 
Environment  (CASTLE).  CASTLE  is  a  modular  shell 
structure  designed  for  simulation  development, 
execution,  and  analysis.  The  modular  design  allows 
CASTLE  to  support  a  wde  range  of  rotary  and  fixed 
wing  aircraft-  (TASTLE  provides  a  standard  looping 
structure,  equations  of  motion,  and  atmospheric  models, 
as  well  as  engineering  analysis  facilities  for  simulation 
validation.^ 

The  data  presented  in  this  paper  is  part  of  an  ongoing 
full  envelope  validation  of  the  V-22  math  model  against 
flight  test  data.  Initially,  time  history  comparisons 
were  used  to  validate  the  dynamic  response  of  the  V-22 
simulation  in  airplane  mode.  For  time  history 


comparison,  the  simulation  inputs  are  overdriven  with 
flight  test  data,  and  the  simulation  outputs  are  plotted 
with  measured  flight  test  outputs.  If  the  simulauon 
outputs  follow  the  flight  test  data  within  a  desired 
tolerance  then  the  simulation  is  validated.  Overall  the 
nonlinear  simulation  matched  the  airplane  mode  flight 
test  data  vtry  well  with  one  exception. 

The  simulation  did  not  match  the  flight  test  data  for 
lateral -directional  maneuvers  with  the  flaps  set  at  twenty 
degrees.  For  this  case,  the  time  history  plots  alone  did 
not  provide  enough  information  to  understand  the 
shortcomings  of  the  simulation.  The  differences 
between  the  simulation  model  and  flight  test  data  for 
these  lateral-directional  maneuvers  were  investigated 
using  linear  model  comparison. 

Linear  model  comparison  is  a  validation  technique  that 
provides  d>Tiamic  mode  and  stability  derivative 
information  which  can  be  used  to  update  the  full 
nonlinear  simulation.^  In  this  technique,  linear  model 
approximations  of  the  simulation  and  airaaft  dynamics 
are  calculated  independently  and  then  compared.  If  the 
simulation  model  matches  the  flight  test  model  within  a 
desired  tolerance  then  the  simulation  is  validated. 

A  number  of  techniques  are  available  to  extract  a  linear 
model  from  a  nonlinear  simulation.^  For  this  analysis', 
the  simulation  linear  model  was  extracted  using  a 
perturbation  technique.  A  second  linear  model 
representing  the  airaaft  dynamics  was  fit  to  flight  lest 
data  using  time  domain  parameter  identification  (PID) 
tools. 

Overview  of  the  V-22  Osprev 

The  V-22  Osprey  is  a  tilt-rotor  aircraft  capable  of  flight 
from  hover  to  high  speed  airplane  mode.  Control  of  the 
V-22  is  accomplished  through  both  conventional 
airplane  and  helicopter  controls,  as  illustrated  by  Figure 
1.  The  airplane  controls  include  an  elevator,  four 
flaperons,  and  two  rudders.  The  helicopter  controls 
include  collective  pilch,  longitudinal  cyclic,  and  lataal 
cyclic  for  each  rotor. 

The  V-22  has  a  digital  fly-by-wire  control  system 
consisting  of  a  primary  flight  control  system  (PFCS) 
and  an  automatic  flight  control  system  (AFCS).  The 
PFCS  provides  pilot  control  input  shaping  and  essential 
feedback  loops.  The  AFCS  is  designed  to  provide  Level 
1  handling  cjualities  and  auto-pilot  functions. 

In  airplane  mode,  the  airaaft  is  primarily  controlled 
with  the  aerodynamic  control  surfaces;  however,  the 
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Roll 


Differential  Collective  Flaperons 

Pitch 

Pitch 


Longitudinal  Cyclic  Elevator 

Yaw 


It  is  common  practice  to  model  the  lateral-directional 
dynamics  of  an  aircraft  with  a  fourth  wder  system/*  ^ 
Typically  in  this  type  of  model,  the  four  stales  are  the 
sideward  velocity,  roll  rale,  yaw  rate,  and  roll  angle. 
The  inputs  to  the  system  are  the  rudder  and  aileron 
deflections. 

To  develop  a  linear  model  based  on  these  states  and 
inputs,  the  equations  of  motion  for  a  rigid  body  are 
significantly  simplified.  The  small  perturbation 
^proach  is  used  to  linearize  the  equations,  and  the 
higher  order  terms  are  neglected.  It  is  assumed  that  the 
lateral-directional  and  longitudinal  equations  of  motion 
may  be  separated."*  Gyroscopic  contributions  from  the 
rotors  are  neglected  because  the  left  and  right  rotors 
rotate  in  opposite  directions.  The  lateral-directional 
dynamic  equations  are  shown  below. 

Fq  +  AF  +  mg^cosGo  =  m(v  -  W^p  +  UQr)  (4) 

A)  ^  ^xxP 
A^q  +  =  ~IxzP 

^  =  /?  +  rtan0Q  (7) 


Differential  Rudders 

Longitudinal  Cyclic 


Thrust 


Collective  Pitch  Throttle 

Command  Command 


For  unaccelerated  initial  conditions,  the  initial  forces 
and  moments  (Yq,  Lq,  and  Nq)  are  zero.  The 
perturbation  force  and  moments  can  each  be 
approximated  with  a  first  order  Taylor  expansion.  The 
partial  derivatives,  Y*,  L.,  and  N*,  are  known  as  the 
dimensional  stability  and  control  derivatives. 

AF  =  n  V  +  Y^p  +  Yj  +  YJ,  +  YJ^  (8) 

AL  =  L^.v  +  LpP  +  Lj  +  Lg^S^  +  Lg^5„  (9) 

AN  =  N^'  +  Npp  +  N^r  +  N^5,+Nj^  (10) 


Figure  1.  V-22  Controls:  Helicopter  (Left)  and 
Airplane  (Right). 


rotor  controls  are  still  significant.  The  PFCS  uses  the 
rotor  cyclic  inputs  to  minimize  the  flapping  angles  of 
each  rotor  and  the  collective  inputs  to  maintain  a 
constant  rotor  speed  and  enhance  thrust  control. 


State-Space  Model  Definition 


The  models  developed  in  this  study  utilize  a  standard 
continuous  time,  state-space  model. 


i  =  Ax  -1-  Bu 
y  =  Cx  +  £hi 

x(^o)  =  Xo 


(1) 

(2) 

(3) 


The"  cross-product  of  inertia,  Ixz,  for  the  V-22  is 
generally  less  than  1%  of  the  moments  of  inertia,  Ixx 
and  Izz.  so  the  roll-yaw  coupling  terms  in  equations  5 
and  6  can  be  ignored  without  significantly  affecting  the 
model  dynamics.  This  yields  the  final  state  equations  for 
the  rigid  body  dynamics. 

m  \m  J  \m  J 

^  ^  (11) 

^COS0O^+— <5,  +— 5, 
m  m 

P  =  ^[L^v  +  LpP+Lrr  +  Lj^  +  Lj^)  (12) 
*xx 
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f  =  -^(N,v+N^p  +  N,r+N^S,  +  NJ.)  (13) 

^72 

0z=p  +  rtan0^  (14) 


fixed  in  airplane  mode,  true  velocity  of  183  knots,  and 
flaps  set  at  twenty  degrees.  A  linear  model  was 
obtained  from  the  full  nonlinear  simulation  using  the 
Linear  Model  Extraction  (LME)  facility.^ 


For  the  V-22,  additional  states  and  inputs  are  required  to 
model  the  dynamics  of  the  rotor  system.  The  slates  are 
the  longitudinal  and  lateral  flapping  angles,  and  the 
inputs  are  the  longitudinal  and  lateral  cyclic  pitch  angles 
and  the  collective  pilch  angles. 

Combining  the  rotor  and  rigid  body  dynamic  models 
generates  an  8th  order  state-space  system.  The  A  matrix 
can  be  divided  into  four  quadrants  as  shown  below, 
where  the  diagonal  quadrants  contain  the  pure  rigid  body 
and  rotor  dynamics,  and  the  off  diagonal  quadrants 
represent  the  cross-coupling  of  the  two  systems. 


A 

^  rigid  body 

A 

^rotor  on  rigid  body 

A 

^ rig  id  body  on  rotor 

A 

^rotor 

The  quadrant  is  defined  by  equations  11-14. 


Similarly,  the  control  matrix,  B,  can  be  separated  into 
quadrants  for  the  control  surface  and  rotor  input  eflects 
on  the  rigid  body  and  rotor  states. 


B 


total 


’  n 

^surfaces  on  rigid  body 

T> 

^ rotor  inputs  on  rigid  body 

D 

^surfaces  on  rotor 

jy 

^rotor  inputs  on  rotor 

(16) 


The  upper  left  quadrant  is  defined  by  equations  11-14. 
Since  the  rigid  body  control  surfaces  have  no  effect  on 
the  rotor  flapping,  the  elements  of  the  lower  left 
quadrant  are  zeros.  The  cyclic  inputs  have  very  little 
effect  on  the  rigid  body  motion,  but  the  collective 
inputs  for  each  rotor  have  a  significant  effect  on  the  yaw 
rale. 

The  C  matrix  is  nearly  identity  since  the  output  vector 
is  approximately  equivalent  to  the  state  vector.  The 
matrix  does  contain  terms  to  transform  the  flapping 
angles  to  the  axes  used  in  flight  test  measurements. 
The  D  matrix  is  identically  zero  since  the  inputs  do  not 
directly  affect  the  outputs. 

Linear  Model  Extraction  from  the  Nonlinear 
Simulation 

The  V-22  simulation  in  the  CASTLE  architecture  was 
initialized  at  the  flight  test  initial  conditions.  The 
initial  conditions  were  level  steady-state  flight,  nacelles 


LME  uses  the  offset  derivative  method  to  extract  a 
linear  model.  In  this  method,  perturbations  are  added  to 
each  input  and  state,  and  the  resulting  changes  in 
outputs  and  state  derivatives  are  recorded.  Integration  is 
frozen  within  the  simulation  so  that  the  state  derivatives 
are  not  allowed  to  propagate.  The  state-space  matrices 
are  then  computed. 


A,= 


Ax,. 


Aup 


”  AX; 


D  = 


(17) 


Aw. 


In  the  above  equation,  i  is  the  stale  derivative  index, 
j  the  state  index,  k  the  control  input  index,  and  n  the 
output  index. 


Since  the  simulation  in  general  is  nonlinear,  the  change 
in  outputs  and  state  derivatives  for  different  size 
perturbations  may  not  be  perfectly  linear.  The  user 
specifies  four  perturbation  sizes  for  each  input  and  stale. 
LME  perturbs  the  simulation  positively  and  negatively 
for  each  of  the  perturbation  sizes  and  picks  the 
perturbation  size  that  results  in  the  best  linearity.  A 
typical  example  of  an  LME  result  for  a  linear 
relationship  is  shown  in  Figure  2. 


Figure  2.  LME  result  for  a  linear  relationship 
between  the  slate  derivative  and 
perturbation  parameter. 
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In  some  cases  the  LME  result  shows  that  a  linear 
relationship  between  an  output  and  a  perturbation  does 
not  exist.  An  example  of  a  nonlinear  relationship  is 
shown  in  Figure  3.  This  figure  shows  that  the  linear 
approximation  of  the  nonlinear  data  is  adequate  near  the 
reference  condition.  The  linearity  information  is  useful 
in  determining  the  quality  of  the  linear  model  with 
respect  to  each  coefficient.  The  output  of  LME  is  the 
state-space  model  and  the  linearity  information  for  each 
coefficient. 


Figure  3,  LME  result  for  a  nonlinear  relationship 
between  the  state  derivative  and 
perturbation  parameter. 


Verification  of  the  LME  Model 

The  state-space  LME  model  was  compared  to  the  full 
nonlinear  simulation  to  examine  the  quality  of  the 
extracted  linear  model.  Two  input  maneuvers,  a  lateral 
stick  ramp  and  a  rudder  doublet,  were  used  to  verify  the 
linear  model.  Each  model  was  driven  with  the  same 
inputs,  and  the  model  outputs  were  compared.  The 
control  inputs  and  the  model  outputs  are  shown  in 
Figure  4  as  a  concatenated  data  set.  The  figure  shows 
that  the  linear  model  accurately  represents  the  lateral- 
directional  dynamics  of  the  full  nonlinear  simulation 
about  the  reference  flight  condition  for  the  two 
maneuvers. 

Linear  Model  Identification  to  Flight  Test 
Data 

Three  lateral-directional  flight  test  maneuvers,  a  roll 
reversal,  a  yaw  reversal,  and  a  yaw  doublet,  were  used 
for  the  parameter  identification.  Visual  inspection  of 
the  30  degree  roll  reversal  showed  good  excitation  of  the 
roll  mode  and  yaw  due  to  roll  characteristics.  The  yaw 


reversal  and  yaw  doublet  showed  good  excitation  of  the 
dutch  roll  mode  on  the  sideward  velocity  and  yaw  rate 
lime  histories.  The  flight  test  data  was  preprocessed  and 
analyzed  for  kinematic  consistency.  The  final  data  set 
showed  good  consistency  with  the  rigid  body  kinematic 
equations.  The  three  maneuvers  were  concatenated  for 
parameter  identification. 

Parameter  identification  was  limited  to  stability  and 
control  derivatives  from  the  rigid  body  lateral-directional 
equations  contained  in  the  upper  left  hand  quadrants  of 
the  A  and  B  matrices.  No  parameters  in  the  rotor 
dynamic  equations  or  the  cross-coupling  matrices  were 
included  in  the  identification.  Identification  of  rotor 
dynamic  parametas  from  maneuvers  with  mainly 
conventional  airplane  control  inputs  was  not  attempted. 
The  LME  state-space  matrices  were  used  as  the  initial 
parameter  values  in  the  identification. 

The  parameters  available  for  identification  formed  a  set 
of  nine  stability  derivatives  and  six  control  derivatives. 
Cramer-Rao  bounds,  which  are  estimates  of  the  standard 
deviation  of  the  identified  parameters,  were  used  to 
assess  the  identifiability  of  each  parameter,^  Three 
stability  derivatives  and  one  control  derivative  had 
relatively  high  Cramer-Rao  bounds  and  were  not 
identified.  These  parameters  were  held  constant  at  the 
LME  model  value. 

The  cost  function  was  defined  as  a  weighted  sum  of 
squared  errors  between  the  measured  outputs,  y ,  and  the 
estimated  outputs,  y .  The  weighting  matrix,  G,  is  a 
diagonal  mauix  that  contains  the  relative  weighting  for 
each  output 

■/“ilclji-y,)’]  (18) 

1=1 

The  rigid  body  states  were  weighted  ten  times  greater 
than  the  rotor  flapping  states  for  the  identification.  The 
rotor  flapping  outputs  were  included  in  the  cost  function 
because  the  rigid  body  equations  and  rotor  dynamic 
equations  are  coupled. 

The  cost  function  was  minimized  through  parameter 
identification  using  the  MATLAB®  computing 
environment.  The  state-space  model  was  overdriven 
with  flight  test  inputs  to  produce  the  estimated  outputs. 
A  Levenburg-Marquardt  routine,  which  is  a  combination 
of  the  Gauss-Newton  and  steepest  descent  methods,  was 
used  to  minimize  the  cost  function.’ 
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Figure  4.  Verification  of  the  LME  model  to  the  full  nonlinear  simulation. 


The  measured  inputs  and  outputs  from  flight  test,  the 
LME  model  outputs,  and  the  PID  model  outputs  ate 
shown  in  Figure  5.  A  marked  improvement  in  the 
identified  model  is  evident  from  the  plots  of  the 
outputs. 

The  improvement  of  the  identification  was 
quantitatively  assessed  using  goodness  of  fit  statistics 
applied  to  the  measured  and  modeled  outputs.  Theifs 
inequality  coefficient  is  a  statistic  that  scales  the  root 


mean  square  error  between  zero  and  one.  An  inequality 
coefficient  value  of  zero  indicates  a  perfect  fit,  and  a 
value  of  one  indicates  the  fit  is  as  bad  as  possible.* 


\ 
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Figure  5.  Hight  test  data  and  LME  and  PID  model  outputs  for  the  parameter  identification  maneuvers. 


Goodness  of  fit  statistics  showing  the  improvement 
from  the  initial  model  to  the  identified  model  are 
presented  in  the  first  part  of  Table  1. 

Verification  of  the  PID  Model 

The  PID  state-space  model  was  verified  against 
indqiendent  flight  test  data  to  assess  the  quality  of  the 
model.  Often  a  model  will  match  the  data  used  in  the 
identification  but  will  not  accurately  predia  independent 
data.  Good  prediction  capability  provides  confidence 


that  the  identified  model  represents  the  actual  aircraft 
dynamics. 

A  yaw  doublet  not  included  in  the  identification  was 
used  to  verify  the  identified  model.  The  identified  model 
was  overdriven  with  the  inputs  from  the  yaw  doublet 
and  the  outputs  were  compared  to  the  flight  test  data. 
Figure  6  shows  the  flight  test  and  model  ouqmts. 
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Table  1.  Goodness  of  fit  statistics  comparing  the 

LME  outputs  and  the  PID  model  outputs  to 
flight  test  data  used  in  the  identification  and 
independent  verification  data. 


U-Theil 

State  Space 
Model 

LME  Model 

PID  Model 

PID  Model 

Outputs 

Identification 
Data  Set 

Identification 
Data  Set 

Verification 
Data  Set 

V 

0.679 

0.079 

0.210 

p 

0.742 

0.072 

0.256 

r 

0.586 

0.079 

0.227 

<t> 

0.491 

0.094 

0.293 

^IR 

0.674 

0.388 

0.411 

^IL 

0.685 

0.310 

0.345 

b,L 

0.728 

0.223 

0.217 

blR 

0.723 

0.145 

0.255 

The  PID  model  adequately  represents  the  verification 
data  at  the  beginning  of  the  maneuver  while  the  rudder 
doublet  is  active.  TTie  model  is  somewhat  overdamped 
during  the  remainder  of  the  maneuver.  Goodness  of  fit 
statistics  for  the  verification  yaw  doublet  are  presented 
in  the  last  column  of  Table  L  As  expected,  the  Theil 
coefficients  for  the  verification  data  set  are  higher  than 
those  for  the  identification  data  set,  but  they  still  show 
an  improvement  over  the  initial  LME  model. 

The  identification  process  produced  a  state-space  model 
with  adequate  dutch  roll  prediction  capability.  Less 
confidence  was  placed  in  the  roll  mode  predictability  for 
two  reasons.  The  thirty  degree  roll  reversal  used  in  the 
identification  stretches  the  linearity  assumption,  and 
indq>endent  flight  test  data  was  not  available  to  verify 
the  roll  response. 

Comparison  of  Linear  Models 

An  obvious  conclusion  that  could  be  drawn  from  Figure 
5  is  that  the  LME  state-space  model  does  not  accurately 
represent  the  dutch  roll  mode  of  the  aircraft.  The  time 
history  plots  for  the  yaw  reversal  from  the  LME  model 
suggest  a  slightly  unstable  dutch  roll  mode. 
Quantifying  this  problem  using  only  a  time  history 
comparison  would  be  very  difficult. 

An  inherent  strength  of  the  linear  model  comparison 
method  is  the  ability  to  compare  the  dynamic  modes 


Table  2,  Lateral-directional  dynamic  stability 
characteristics  for  the  LME  and  PDD 
models. 


Dutch  Roll 
Mode 

Spiral 

Mode 

Roll 

Mode 

% 

Ts 

Tr 

LME 

Model 

-0.004 

0.820 

7.94 

1,26 

PID 

Model 

0.041 

1.005 

8.87 

1.37 

Table  3.  Nondimensional  stability  derivative  values 
from  the  LME  and  PID  models. 


Stability 

Derivative 

LME 

Model 

PE) 

Model 

Diffaence 

Cyp 

-2.41 

-2.42 

-0.01 

Cyp 

-0.464 

-0.464 

No  ID 

Cy, 

0.506 

0.506 

No  ID 

Cy5r 

-0.236 

-0.240 

-0.004 

Cysa 

-0.226 

-0.226 

No  ID 

Cip 

-0.342 

-0.406 

-0.064 

Qp 

-1.20 

-1.18 

0.02 

Cu 

-0.169 

-0.169 

NoE) 

Ci5r 

-0.017 

-0.072 

-0.055 

Ci5a 

-0.350 

-0.320 

0.030 

0.042 

0.138 

0.096 

Cnp 

-0.640 

-0.683 

-0.043 

^Nr 

-1.77 

-1.85 

-0.08 

CnBt 

0.173 

0.159 

-0.014 

Civsa 

-0.0217 

-0.0400 

-0.0183 

contained  in  the  state-space  models.  The  eigenvalues  of 
the  A  matrix  which  describe  the  modes  of  the  aircraft  are 
easily  calculated  in  MATLAB*.  The  eigenvalues  of  the 
lateral-directional  rigid  body  equations  are  shown  on  a 
complex  plane  plot  in  Figure  7.  The  damping  ratio  and 
undamped  natural  frequency  of  the  dutch  roll  mode, 
along  with  the  time  constants  of  the  spiral  and  roll 
modes,  are  shown  in  Table  2.  The  acceptable  tolerance 
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Figure  6.  Right  test  data  and  FED  model  outputs  for  an  independent  verification  yaw  doublet. 

for  error  in  the  mode  comparison  depends  on  the  specific  the  dutch  roll  natural  fiequency  of  the  PE)  model, 

aircraft  and  simulation  application.  A  general  tolerance  Based  on  the  proposed  tolerance,  the  LME  model  is  not 

of  5%  has  been  proposed  for  this  type  of  analysis.^  a  valid  representation  of  the  dutch  roll  mode  exhibited  in 

flight  test. 

Both  the  LME  and  PE)  models  contain  classic  second 

OTdo'  dutch  roll  modes.  The  dutch  roll  mode  of  the  The  LME  and  PE)  models  both  contain  first  order  spiral 

LME  model  is  slightly  unstable  as  evident  from  the  and  roll  modes.  The  spiral  mode  time  constant  of  the 

negative  damping  ratio.  The  dutch  roll  damping  ratio  of  LME  model  is  10%  less  than  that  of  the  PE)  model, 

the  identified  model  is  larger  in  magnitude  and  positive.  The  roll  mode  time  constant  of  the  LME  model  is  8% 

Ihe  LME  dutch  roll  natural  frequency  is  18%  less  than  less  than  the  PE)  model  value. 


9 

American  Institute  of  Aeronautics  and  Astronautics 


1.0 

0.8 

0.6 

0.4 

0.21 


Imag 

Axis 


-0.2 

-0.4 

-0.6 

-0.8 


-1.0 


X  =  LME  moctei 
K  =PIDrTKXlel 


-1.5 


Dutch  Roll 


Roll 


Spiral 


Dutch  Roll 


-1.0 


-0.5 

Rea!  Axis 


0.5 


Figure  7,  Complex  plane  plot  of  the  lateral- 

directional  dynamic  modes  from  the  LME 
and  the  FED  models. 


In  addition  to  the  dynamic  mode  information,  the 
aircraft  stability  and  control  derivatives  can  be  calculated 
from  the  state-space  models.  The  dimensional  stability 
derivatives  are  imbedded  in  the  elements  of  the  state- 
space  model  as  shown  in  equations  11-13.  The  rigid 
body  stability  and  control  derivatives  wae  extracted 
from  the  A  and  B  matrices  and  nondimensionalized. 
Table  3  contains  the  nondimensional  stability  and 
control  derivatives  calculated  from  the  LME  and 
identified  models. 

Simulation  Stability  Derivative 
Investigation 

Two  approaches  can  be  taken  to  incorporate  the 
increments  in  stability  and  control  derivatives  into  the 
nonlinear  simulation.  The  most  straightforward 
approach  is  to  directly  acid  the  increments  into  the 
aerodynamic  tace  and  moment  equations.  This 
approach  is  only  appropriate  for  simple  simulation 
models.  The  nonlinear  nature  of  the  V-22  simulation 
does  not  lend  itself  to  this  direct  method.  For  complex 
models  such  as  the  V-22,  the  stability  and  control 
derivative  information  can  be  used  to  locate  weaknesses 
in  the  physically  based  equations  of  the  simulation 
model. 

The  shaded  stability  derivative  in  Table  3,  showed 
the  greatest  difference  between  the  LME  and  FID 


models.  is  the  nondimensional  yawing  moment 
coefficient  due  to  sideslip.  must  be  positive  for  the 
aircraft  to  exhibit  directional  stability.  In  the  nonlinear 
simulation,  this  coefficient  contains  contributions  from 
the  aerodynamic  and  rotor  models.  The  signs  and 
magnitudes  of  these  contributions  can  be  examined  in 
the  LME  environment. 

Examination  of  the  components  showed  that  the 
rotor  model  contribution  is  destabilizing  (negative)  and 
the  aCTodynamic  model  contribution  is  stabilizing 
(positive).  The  rotor  model  contribution  is  only 
slightly  smaller  in  magnitude  than  the  aerodynamic 
model  contribution.  The  rotor  and  aerodynamic  models 
combined  to  produce  a  total  simulation  that 
provides  weak  directional  stability. 

Investigation  into  the  rotor  model  equations  showed  that 
the  magnitude  of  the  rotor  model  contribution  is 
dependent  on  a  sideward  velocity  effect  model.  This 
sideward  velocity  effect  is  intended  to  model  the 
interference  effect  that  one  rotor  would  have  on  the 
other.  For  high  speed  airplane  mode  flight,  physical 
intuition  suggests  that  the  interference  effect  between 
the  rotors  should  be  negligible.  Elimination  of  the 
sideward  velocity  effect  in  the  nonlinear  simulation 
decreased  the  rotor  destabilizing  effect  on  the  directional 
stability,  and  the  total  moved  closer  to  the  FID 
predicted  value. 

Further  investigation  is  required  to  determine  the  history 
of  the  sideward  velocity  effect  and  assess  the  validity  of 
the  effect  for  airplane  mode  flight.  In  this  investigation 
the  linear  model  comparison  technique  has  provided 
valuable  insight  into  the  validation  process  and  may 
have  located  a  weakness  in  the  rotor  model  physical 
equations. 


Conclusions 

A  linear  model  was  extracted  from  the  V-22  full 
nonlinear  simulation  using  LME.  This  model  was 
verified  against  time  history  data  from  the  nonlinear 
simulation.  The  linear  model  accurately  represented  the 
lateral-directional  dynamics  of  the  verification 
maneuvers.  A  second  linear  model  was  fit  to  flight  test 
data  using  parameter  identification  tools.  This  model 
was  verified  using  an  independent  maneuver  not  included 
in  the  identification  process.  The  linear  model  exhibited 
adequate  predictability  of  the  dutch  roll  mode  dynamics. 
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sufficient  data  for  identification  and  verification. 
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